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PINCHED  PROPAGATION  OF  HIGH-POWER, 

PULSED  ELECTRON  BEAMS  FOR  WELDING  AND  MATERIALS 

PROCESSING  APPLICATIONS 

I.  INTRODUCTION 

Many  industrial  processes  require  that  a  localized  heat  source  be  applied  to  a  metal, 
ceramic,  or  other  material.  These  processes  include  welding,  bonding,  cutting,  annealing,  surface 
hardening,  and  fabrication  of  advanced  materials.  Although  conventional  heat  sources  are 
commonly  employed  for  these  purposes,  lasers  and  electron  beams  have  also  been  used. 
Electron  beams  are  especially  promising  because  they  offer  high  irradiance,  the  ability  to  pass 
through  metal  vapors  with  little  degradation,  and  high  "wallplug-to- workpiece"  efficiency.1  For 
example,  compared  with  traditional  arc  welders,  electron  beams1'9  can  produce  narrower  welds  at 
much  higher  speeds  and  lower  heat  input.1'3 

Nevertheless,  electron  beam  processing  is  not  widely  used  commercially,  except  in  the 
automotive  industry  in  the  United  States  where  atmospheric  electron  beam  welding  has  been 
used  since  the  mid-sixties.3  In  part,  this  is  due  to  the  limitations  of  the  technology  that  has  been 
employed  to  date.  For  example,  present-day  commercial  electron  beam  welders  typically  operate 
at  energies  less  than  200  keV  and  at  average  currents  less  than  1  amp.  When  operating  in  the 
atmosphere  at  these  energies,  scattering  of  the  electrons  off  air  molecules  causes  rapid  expansion 
of  the  beam.  To  reduce  beam  divergence,  magnetic  focusing1 3,4  is  deployed  at  the  point  where 
the  beam  is  released  into  the  atmosphere.  Beam  expansion  can  be  further  reduced  by  taking 
advantage  of  the  reduction  of  air  density  within  the  beam-heated  air  channel  ("hole-boring"),5,9 
or  by  propagating  the  beam  within  a  helium  jet  rather  than  open  air.2,5  Nevertheless,  the  depth  of 
focus  remains  limited,  and  the  beam  source  must  be  operated  within  a  few  cm  of  the  workpiece. 
A  larger  standoff  would  be  far  more  attractive  because  it  would  ease  target  setup,  produce  deep 
and  parallel  weld  joints,  and  minimize  variations  in  beam  focus  that  result  from  power-supply 
ripple,  fluctuations  in  the  weld-cavity  depth,  or  curvature  of  the  workpiece  surface.4,6  Perhaps 
most  significantly,  splash-back  of  debris  into  the  accelerator  would  be  reduced;  at  present, 
splash-back  is  a  major  problem6  due  to  arcing  within  the  accelerator.  To  minimize  arcing4  (and 
also  to  minimize  the  need  for  x-ray  shielding5),  beam  energies  greater  than  200  keV  have 
generally  been  avoided. 

Many  of  these  limitations  can  be  circumvented  by  operating  in  vacuum,  but  this 

profoundly  changes  the  processing  conditions  and  is  unacceptable  for  certain  applications.  At 
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best,  the  use  of  a  vacuum  system  increases  the  cost  considerably,  both  in  terms  of  equipment 
needed  and  operating  time  lost  to  chamber  evacuation.6 

Another  consideration  is  that  at  these  energies,  the  electrons  penetrate  only  a  short  distance 
into  the  workpiece  (10-100  microns  into  steel,  for  example),  so  that  the  beam  deposition  is 
confined  to  the  surface.  As  is  the  case  with  laser  welding  and  conventional  heat  sources,  the 
interior  is  heated  through  a  combination  of  thermal  conduction  and  convection.  However,  these 
processes  reduce  system  efficiency  by  spreading  the  energy  to  areas  outside  the  target  zone,  and 
they  produce  a  "nail-head"  profile  in  which  the  weld  radius  decreases  with  depth.7  Furthermore, 
the  weld  joint  is  vulnerable  to  irregularities  in  the  target  material  and  the  cleanliness  of  its 
surface.6,8  An  additional  problem  specific  to  electron  beams  is  joint  detracking,  caused  by 
magnetic  beam  deflection  due  to  magnetic  materials,  asymmetrical  currents  in  the  liquid  weld 
pool,  or  thermoelectric  currents  produced  by  the  heating  of  dissimilar  metals.4,6,8 

As  has  been  recently  pointed  out,10  many  of  these  problems  can  be  eliminated  by  using 
high  current,  pulsed  electron  beams  with  energy  in  the  relativistic  range  ( >  1  MeV).  For  the  sake 
of  brevity,  we  shall  adopt  the  terminology  of  a  recent  workshop11  on  electron-beam  processing, 
and  refer  to  beams  of  this  type  as  high-energy  electron  beams  (HEEB).  Pulsed  rather  than  dc 
operation  allows  much  higher  accelerating  voltages  to  be  applied  without  insulator  breakdown,12 
and  research  laboratories  routinely  utilize  this  technology  to  produce  electron  beams  with 
energies  of  several  MeV,  currents  of  1  kA  or  more,  and  pulse  lengths  of  10-1000  ns.  The  beams 
have  been  generated  using  diodes,  betatrons,  and  staged  induction  accelerators.13'18  Reliable  and 
reproducible  operation  has  been  achieved  for  many  millions  of  shots  and  at  pulse  repetition  rates 
up  to  1  kHz.16  Although  x-ray  production  is  much  higher  at  these  energies,  adequate  shielding 
has  not  been  a  problem.  In  fact,  the  production  of  copious  x-rays  offers  the  potential  for  real¬ 
time,  in-situ  diagnostics.  However,  for  typical  materials  processing  applications,  the  beam 
energy  should  be  restricted  to  less  than  10  MeV,  in  order  to  avoid  nuclear  activation  of  long-lived 
radioactive  isotopes  in  the  target.  It  should  be  noted  that  low-current,  non-pulsed  electron  beams 
in  this  energy  range  are  already  in  use  for  certain  industrial  materials  processes.19'21 

Pulsed  beams  with  high  current  and  high  electron  energy  offer  significant  advantages  over 
the  dc  low-power  beams  used  presently.  The  first  class  of  advantages  relates  to  propagation  to 
the  workpiece.  High  electron  energy,  in  itself,  makes  the  beam  "stiff'  and  thus  less  sensitive  to 
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stray  magnetic  fields  or  gas  scattering.22  This  reduces  joint  mistracking  and  permits  a  somewhat 
increased  standoff  distance.  More  importantly,  beams  with  both  high  energy  and  high  current  (in 
the  kA  range)  can  propagate  in  dense  gas  in  a  self-pinched  mode,  because  the  beam  space  charge 
is  neutralized  by  beam-ionized  plasma.  As  we  shall  show,  self-pinching  permits  stable,  tight 
propagation  through  tens  of  cm  or  more  of  ambient  air.  Adequate  standoff  can  then  be  provided, 
without  need  of  vacuum  pumping.  A  second  class  of  advantages  relates  to  the  very  different  way 
in  which  a  HEEB  interacts  with  the  workpiece.  High-energy  beams  deposit  energy  nearly 
instantaneously  and  deep  into  the  target  material.  For  example,  the  range  of  1-10  MeV  electrons 
in  dense  materials  like  steel  or  titanium  is  approximately  0.1-1  cm.  This  deep  penetration  makes 
material  processing  much  less  dependent  on  the  condition  and  cleanliness  of  the  surface. 
Furthermore,  because  the  energy  is  deposited  in  a  time  much  shorter  than  the  hydrodynamic  or 
thermal  conduction  times  on  which  the  material  responds,  the  heating  is  more  efficient  and  the 
spatial  extent  of  the  heat-affected  zone8  is  reduced,  which  can  be  important  for  applications  like 
welding.  For  applications  like  hardening,  it  may  be  useful  to  use  pulse  energy  densities  that  are 
large  enough  to  drive  localized  shocks  as  well  as  heat  the  material.  For  other  applications,  it  may 
be  convenient  to  vary  the  beam  parameters  during  the  pulse  train,  to  tune  the  time  dependence  of 
the  energy  deposition  (as  needed,  for  example,  when  making  continuous  weld  joints4).  In  short, 
high-energy,  high-current,  pulsed  electron  beams  offer  new  capabilities  and  the  potential  to 
overcome  many  of  the  limitations  of  the  commercial  devices  in  use  today. 

A  large  number  of  applications  for  HEEB  welding  and  materials  processing  has  been 
suggested,"  but  the  beam  parameters  required  for  these  applications  are  diverse  and  not  well 
known  at  the  present  time.  Some  applications,  such  as  welding,  require  small  beam  radii 
(ab  ~  0.1  cm)  and  high  source  brightness  from  the  accelerator.  Other  applications,  such  as 
bonding,  surface  modification,  and  materials  fabrication,  typically  require  beams  with  larger  radii 
and  lower  brightness,  but  possibly  with  larger  total  current.  This  paper  addresses  the  limits 
imposed  on  the  choice  of  HEEB  parameters  by  the  physics  of  beam  propagation  in  the 
atmosphere  or  other  dense  gas,  over  distances  of  a  few  tens  of  cm  from  the  accelerator  to  the 
target.  For  specificity,  we  further  restrict  the  analysis  to  a  "welding"  regime  where  aj,  ~  0.1  cm, 
and  a  "materials  processing"  regime  where  we  arbitrarily  assume  that  at,  ~  1  cm.  Everything  else 
being  equal,  it  is  always  easier  to  propagate  a  beam  with  a  larger  radius. 
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The  propagation  of  high-power  electron  beams  is  a  complic  'd  field  that  has  been  the 
subject  of  an  extensive  literature.  To  analyze  the  various  physics  issues,  we  draw  on  much  of 
this  previous  work.  However,  we  shall  point  out  that  there  are  significant  physical  differences 
between  the  parameter  regime  of  interest  here  (e.g.,  for  welding,  beam  currents  Ib  -  1  kA,  beam 
radius  ab  ~  0. 1  cm,  propagation  distance  zt  ~  30  cm)  and  the  focus  of  most  of  the  earlier  work 
(Ib>  10  kA,  larger  ab,  and  very  long  zt).  In  Sec.  II,  we  review  HEEB  generation  and  the  physics 
of  self-pinching  in  dense  gas.  In  Sec.  Ill  we  discuss  the  resistive  hose  instability  and  its  impact 
on  beam  propagation  and  energy  deposition  at  the  target;  several  stabilizing  effects,  including 
phase-mixing,  gas  scattering,  and  pipe  guiding,  are  shown  to  be  important.  In  Sec.  IV,  we 
illustrate  the  propagation  limitations  on  HEEB  using  both  analytic  examples  and  numerical 
simulations  computed  using  a  three-dimensional,  fully  electromagnetic  particle  code  (SARLAC) 
developed  for  studying  the  propagation  of  relativistic  electron  beams  in  the  atmosphere. 
Miscellaneous  issues  such  as  gas  heating,  deflection  from  stray  fields,  and  accelerator  interface 
problems  are  discussed  in  Sec.  V.  The  detailed  mathematical  treatment  of  hose  instability  is 
considered  in  the  Appendix. 
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II.  REVIEW  OF  ELECTRON  BEAM  PROPAGATION  PHYSICS 


A.  Beam  Production  and  Conditioning 

The  production  of  high-power  electron  beams  utilizes  pulsed  power  technology  that  has 
been  under  development  since  the  1960’s.  In  the  most  common  configuration,  a  prime  power 
supply  is  connected  to  a  Marx  generator  consisting  of  a  series  of  capacitors  that  are  charged  in 
parallel  and  discharged  in  series.  The  output  from  the  Marx  generator  is  usually  passed  to  an 
intermediate  storage  capacitor,  followed  by  a  pulse  forming  line  (PFL).  The  high  voltage,  short- 
pulse  output  from  the  PFL  is  then  switched  into  a  load.  If  the  load  is  a  diode,  electrons  are 
produced  at  the  cathode  and  accelerated  towards  a  positive  anode.  High-power  diodes  are 
generally  single-stage  devices  inherently  limited  to  a  few  megavolts  or  less,  and  in  the  past  they 
have  usually  been  single-pulse  research  machines  incapable  of  the  high  throughput  needed  for 
commercial  applications.  However,  there  has  been  recent  interest  in  developing  high  average 
power,  pulsed  diodes  for  applications  like  high-power  microwaves.23 

Higher  voltages  can  be  generated  by  deploying  a  series  of  linear  induction  accelerating 
cavities.  This  is  the  basis  of  several  large  high-current  accelerators,  including  the  RADLAC 
accelerator  at  Sandia  National  Laboratories13  and  the  Experimental  Test  Accelerators  (ETA  and 
ETA-II)  and  the  Advanced  Test  Accelerator  (ATA)  at  Lawrence  Livermore  National 
Laboratory.16'18  ETA-II  is  particularly  of  interest  because  its  brightness,  pulse  repetition  rate, 
current,  and  energy  are  in  an  appropriate  regime  for  welding  and  materials  processing.  However, 
all  of  these  devices  are  large  and  expensive,  and  their  limited  number  of  pulses  in  a  burst  makes 
the  average  power  insufficient  for  most  commercial  applications. 

New  accelerator  concepts  are  currently  being  developed  that  may  soon  be  capable  of 
delivering  high  average  power  at  a  cost  low  enough  to  be  of  commercial  interest.  The  SNOMAD 
device,12  which  is  being  developed  at  Science  Research  Laboratory,  uses  solid-state  switches  to 
produce  0.5  kA,  0.5  MeV  beams  with  a  pulse  repetition  frequency  (PRF)  up  to  5  kHz. 
Additional  modules  will  raise  the  energy  to  1.5  or  2.5  MeV  in  the  near  future.  The  Repetitive 
High  Energy  Pulsed  Power  (RHEPP)  system12  being  developed  at  Sandia  National  Laboratories, 
uses  a  magnetic  pulse  compressor  to  drive  a  linear  induction  voltage  adder  with  a  diode  load. 
This  system  is  designed  to  accelerate  either  electron  or  ion  beams,  with  higher  currents  (~25  kA) 
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but  lower  PRF  (~120  Hz)  than  SNOMAD.  The  pulsed  power  system  operates  at  0.9  MV,  with 
plans  to  produce  a  2.5  MeV  beam. 

Previous  studies  of  long-range  beam  propagation  have  focused  on  the  development  and 
use  of  post-accelerator  conditioning  cells24-26  to  suppress  the  resistive  hose  instability.  We  shall 
discuss  this  beam-plasma  instability  in  Sec.  IE,  and  shall  show  that  it  is  a  matter  of  some  concern 
for  the  present  applications.  However,  we  shall  conclude  that  the  use  of  beam  conditioning  is  not 
necessary  or  even  useful  for  the  relatively  short-range  propagation  considered  here. 

B.  Seif-Pinched  Propagation 

In  using  HEEB  for  welding  and  materials  processing,  we  wish  to  deliver  a  beam  of  high 
current  density  onto  a  target  piece  located  several  tens  of  cm  from  the  accelerator  in  ambient  air 
or  other  dense  gas.  On  its  way  to  the  target,  the  beam  collides  with  gas  molecules,  expands  as  a 
result  of  its  transverse  pressure,  and  may  suffer  electromagnetic  instabilities.  At  some  point, 
these  processes  degrade  the  beam  to  the  extent  that  it  is  no  longer  useful.  In  this  and  the 
following  sections,  we  review  the  physics  of  beam  propagation  in  dense  gas  and  determine  the 
constraints  that  propagation  imposes  on  HEEB  applications. 

We  simplify  the  analysis  by  assuming  the  beam  to  be  highly  relativistic  and  paraxial, 
meaning  that  the  beam  electrons  travel  in  the  z  direction  with  an  axial  velocity  vz  «  c  and  a 
transverse  velocity  \±  «  vz.  These  assumptions  decouple  the  transverse  motion  of  the  beam 
from  its  longitudinal  motion,  and  permit  us  to  describe  beam  expansion  and  contraction  in  terms 
of  a  transverse  beam  temperature,  T±=  <ymv±/2>,  where  m  is  the  electron  rest  mass  and 
7  =  (l-v2/c2)',/2  »  1  is  the  relativistic  mass  factor.  In  what  follows,  we  shall  assume  that  7  is 
the  same  (mono-energetic)  for  all  electrons  in  a  given  beam  slice.  These  assumptions  are 
appropriate  for  HEEB,  but  not  for  the  low-energy  beams  conventionally  used  for  welding  and 
materials  processing. 

All  beams  are  produced  with  a  finite  temperature  Tx,  and  this  temperature  increases  as  the 
beam  collides  with  gas  molecules  on  its  way  to  the  target.  In  the  absence  of  a  radial  pinch  force, 
the  dispersion  of  velocity  angles  associated  with  T±  causes  the  beam  to  expand  until  it  is 
ultimately  too  tenuous  for  the  desired  application.  The  useful  range  can  be  extended,  however, 
by  providing  a  pinch  force.  For  example,  a  solenoidal  magnetic  field  can  be  applied  to  confine 
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the  beam,  although  this  is  often  impractical  because  of  geometrical  constraints  or  the  cost  of  the 
field  coils. 

If  the  beam  is  sufficiently  intense  and  propagating  in  dense  gas,  a  self-pinch  force  develops 
naturally  as  the  beam  ionizes  the  gas  and  creates  a  conducting  plasma.  The  plasma  fully 
neutralizes  the  beam  space  charge  but  is  too  resistive  to  fully  neutralize  the  beam  current.  As  a 
result,  a  finite  net  current  (equal  to  the  sum  of  the  beam  current  and  the  plasma  return  current)  is 
left  flowing  parallel  to  the  beam.  The  net  current  pinches  the  beam  magnetically  and  retards  its 
expansion.  The  plasma  thus  enables  the  beam  to  propagate  in  self-pinched  equilibrium  without 
external  guide  fields.  Self-pinched  propagation  is  the  feature  that  makes  HEEB  welding  and 
materials  processing  possible  with  significant  standoff  in  the  atmosphere.  Figure  1  provides  a 
conceptual  illustration  of  self-pinched  propagation  of  an  intense  beam  to  a  target. 

In  self-pinched  equilibrium,  the  magnetic  pinch  force  from  the  net  current  balances  the 
transverse  beam  pressure: 

h  'Vi1  *  -eVB«  =  -eV  — rc  •  111 


Here  nb  is  the  beam  density,  e  is  the  electron  charge,  p  ■  1  is  the  axial  beam  velocity  normalized 
to  the  speed  of  light  c,  B9  is  the  azimuthal  magnetic  field,  and  in(r)  is  the  net  current  enclosed 
within  radius  r.  Integrating  Eq.  (1)  twice  yields  the  Bennett  pinch  condition  for  the  average  beam 
temperature:27 


epi 


eff 


2c 


where 
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eff 
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n 


(r)  . 


(2) 


(3) 


Here  Ib  is  the  total  beam  current,  ib(r)  is  the  beam  current  within  radius  r,  and  Ieff  is  termed  the 
effective  pinch  current.  If  there  is  no  plasma  return  current,  in(r)  =  ib(r)  and  Ieff = Ib. 


7 


Lee27  has  shown  that  gas  scattering  causes  a  self-pinched  beam  to  evolve  to  an  isothermal 
state,  Tx  =  Tb  for  all  r,  in  which  the  beam  density  assumes  a  Bennett  profile: 
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(r)  = 


„  2  2  2  2 
e/Jcira^  (1+r  /a^) 


(4) 


where  ab  is  the  half-current  radius.  In  this  equilibrium  state,  the  individual  beam  electrons  are 
not  localized  at  fixed  r  but  undergo  sinusoidal-like  betatron  trajectories  about  the  beam  axis,  with 
a  characteristic  axial  wavelength  given  by 
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where  IA  =  07mc3/e  =  1707  kA  is  the  Alfven  current.  The  distance  X^  is  termed  the  average 
betatron  wavelength,  and  it  characterizes  the  propagation  distance  needed  for  equilibration.  For 
paraxial  beams,  Xtf »  a,,. 


C.  Elastic  Gas  Scattering 

Elastic  scattering  of  beam  electrons  off  gas  molecules  eventually  leads  to  beam  expansion, 
even  if  the  beam  is  self-pinched  (but  much  more  rapidly  if  the  beam  is  unpinched).  Scattering 
acts  as  a  source  of  beam  transverse  energy,  given  by22 
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where  na  is  the  atomic  number  density,  Za  is  the  atomic  number  of  the  gas,  and  z  is  propagation 
distance  from  the  beam  source.  If  the  beam  is  unpinched,  its  radius  increases  with  z  as28 
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where  a0  is  the  initial  beam  radius,  0f  is  the  focusing  angle,  0t  =  (2To/027mc2),/2  is  the  initial 
transverse  thermal  spread,  and  T0  is  the  temperature  at  the  injection  point  (z  =  0)  into  the  gas. 
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For  unpinched  beams,  the  radius  can  be  minimized  at  target  location  z,  by  using  a  focusing  angle 
Of  -  a^z,,  as  can  be  shown  by  setting  the  derivative  of  Eq.  (7a)  with  respect  to  Of  equal  to  zero. 
The  radius  at  z,  is  then  given  by 
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Thus  a  minimum  beam  spot  size  ab(zt)  is  set  by  both  the  initial  thermal  spread  and  the  additional 
velocity  spread  resulting  from  gas  scattering.  For  zt  =  30  cm  in  atmospheric  air,  gas  scattering 
alone  restricts  the  radius  to  ab(zl)  >  1  cm  for  7  £  21  (10  MeV). 

In  Fig.  2  we  plot  ab(z)  for  an  unpinched  beam  propagating  in  full  density  air  with 
=  0. 1  cm,  optimum  focusing  angle  0{  =  0.003,  initial  thermal  spread  0t  =  0.002,  and  energy 
e0  =  3,  5  and  10  MeV.  As  is  evident,  the  beam  remains  narrow  for  only  a  few  cm  of 
propagation,  expanding  to  ab  >  1  cm  at  z  2s  30  cm.  While  such  radii  might  be  acceptable  for 
materials  processing,  they  arc  far  in  excess  of  the  typical  requirement  for  welding, 
ab(zt)  “  0.1  cm.  The  bottom  (dashed)  curve  in  Fig.  2  shows  the  propagation  of  the  same  10  MeV 
beam  in  vacuum  (na  =  0).  In  this  case,  the  beam  focuses  to  a^z,)  =  0tz,  =  0.06  cm  at  Zj  =  30  cm. 
Equations  (6)  and  (7)  indicate  that  beam  expansion  in  air  can  be  decreased  by  directing  a  jet  of 
gas  of  low  atomic  number  (e.g.,  helium)  in  the  path  of  the  beam,2,5  or  by  allowing  the  beam  to 
heat  the  air  and  thus  "bore"  a  reduced-density  channel  in  it.9  Both  techniques  have  been  shown 
to  improve  the  performance  of  conventional  electron  beam  welders  operating  in  full  density  air, 
but  the  performance  is  still  inferior  to  that  achieved  in  vacuum.3,6 

The  presence  of  a  pinch  force  reduces  the  expansion  rate  and  changes  the  nature  of  the 
expansion  process.  If  the  beam  power  is  sufficient,  the  beam  maintains  a  self-pinched 
equilibrium  satisfying  Eq.  (2).  As  scattering  adds  to  the  beam’s  transverse  thermal  energy,  the 
beam  expands  adiabatically  to  maintain  the  pressure  balance  condition,  Tx  =  TB.  The  beam 
density  then  decreases  exponentially  according  to 


n,  3z 


(8) 
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The  beam  radius,  which  scales  as  nb  l/2  by  particle  conservation,  therefore  undergoes  "Nordsieck 


expansion 


«27 


a,  (z )  =  a  exp (z/L  ) 

b  o  N 


(9) 


where  the  Nordsieck  length  LN  =  2TB/(dTi/dz)s  is  proportional  to  TlCff t0  lowest  order-  A  more 
careful  analysis  by  Hughes  and  Godfrey29  has  shown  that  large-angle  scattering  weakens  the 
dependence  of  LN  on  7lcff,  and  yields  an  approximate  expression  in  the  HEEB  parameter  regime. 
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where  Ieff  is  in  kA  and  p  =  na/5.4xl019  cm'3  is  the  atomic  number  density  in  atmospheres  (atm)  of 

STP  air.  In  air,  Za  =  7.2  so  that  LN  «  6Ap-\y\K)M  cm. 

Equations  (9)  and  (10)  show  that  7leff  must  be  large  to  limit  beam  expansion  due  to  gas 

scattering.  It  is  convenient  to  define  an  expansion  factor 

g  =  a  /a  ,  (11) 

a  max  o 

where  a^  is  the  maximum  beam  radius  allowed  at  the  target.  The  Nordsieck  length  must  then 
satisfy  Lfi  >  z/ln(ga),  so  that  7lcff  must  exceed  a  minimum  value, 
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As  an  illustration,  consider  a  beam  that  propagates  a  distance  zt  =  30  cm  in  atmospheric  air 
(p  =  1).  If  we  require  ga  =  f7,  then  it  is  necessary  that  LN  ^  3z,  =  90  cm,  and  therefore  that 
7Ieff^  20  kA.  Ieff  cannot  exceed  Ib  for  self-pinched  beams,  and  we  shall  see  in  Sec.  II-E  that 
Ieff « Ib  for  the  cases  of  interest.  Thus  the  pulsed  beam  power  P  must  satisfy  7Ib  *  20  kA,  or 
roughly  P>  10  GW.  If  radiation  safety  limits  the  beam  energy  to  10  MeV  (7  *  21),  the  beam 

current  If,  must  exceed  1  kA. 

Figure  3  shows  plots  of  a^z)  from  Eqs.  (9)  and  (10),  for  a  beam  injected  with  a0  =  0.1  cm, 
=  0,  Tj_  matched  to  a  pinch  current  of  Ieff  =  1  kA,  and  at  energies  of  e0  =  3,  5  and  10  MeV  in 


full  density  air.  The  expansion,  although  often  appreciable,  is  much  slower  than  in  the  unpinched 
case  plotted  in  Fig.  2. 

We  conclude  from  this  discussion  that  high  7  alone  (10  MeV)  can  keep  the  beam  radius  to 
under  a  few  cm  over  propagation  distances  of  a  few  tens  of  cm  in  atmospheric  air  or  like  gases. 
This  may  be  sufficient  for  some  materials  processing  applications.  However,  to  keep  the  beam 
radius  well  under  1  cm  for  applications  like  welding,  self-pinched  propagation  is  needed.  For 
propagation  through  30  cm  in  full  atmosphere,  a  minimum  current  and  power  of  1  kA  and 
10  MW  are  typically  required.  We  shall  use  the  latter  values  as  our  nominal  HEEB  parameters. 

Equations  (9)  and  (10)  for  self-pinched  propagation  govern  only  if  the  characteristic  length 
for  the  beam  to  establish  pressure  balance  (T±  =  TB)  is  shorter  than  the  characteristic 
expansion  length  LN, 

*0  <  ln  •  <13a) 

This  condition  is  equivalent  to  the  statement  that  self-pinched  expansion  occurs  only  if  it  is 
slower  than  free  expansion.  Using  Eqs.  (5)  and  (10),  we  can  rewrite  (13a)  as  a  condition  on  the 
beam  radius. 
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where  Icff  is  again  in  kA.  In  atmospheric  air,  condition  (13b)  gives  ^  <  1  cm  for  Ieff  =  1  kA  and 
7  £  2 1 .  Recall,  however,  that  self-pinching  may  not  be  needed  if  radii  at  the  target  in  excess  of 
1  cm  are  adequate  for  the  application. 


D.  Beam  Energy  Loss  and  Erosion 

Collisions  with  the  air  molecules  not  only  scatter  the  beam  but  also  extract  energy  from  the 
beam.  For  beam  energies  e  £  10  MeV,  the  dominant  loss  comes  from  inelastic  collisions  at  a 
rate22 
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where  W0  is  a  typical  molecular  excitation  energy.  In  atmospheric  air,  if  we  use  the  estimate 
W0«  lOeV, 

^  “  [2  +  0.4  ln(7)]  keV/cm.  (14b) 

az 

The  collisional  energy  loss  over  30  cm  of  propagation  is  thus  ~  0. 1  Me  V,  which  is  negligible  for 
beam  energies  of  several  MeV. 

The  beam  also  loses  energy  to  the  inductive  electric  field  that  arises  in  response  to  the 
beam  current.  Averaged  over  the  beam  pulse,  this  field  satisfies 

E  S  LI  /T  ,  (15) 

z  bp 

where  rp  is  the  pulse  duration  and  L  ~  10  nH/cm  is  the  distributed  circuit  inductance  discussed  in 
Sec.  II-E.  For  rp  *  10  ns  and  I5  =  1  kA,  the  time-averaged  field  is  1  kV/cm  or  less.  The  beam 
energy  lost  over  30  cm  of  propagation  is  therefore  less  than  30  keV  per  electron,  which  is 
negligible.  In  practice,  the  ohmic  losses  are  confined  mainly  to  the  beam  head  rather  than 
uniformly  distributed  over  the  pulse.  This  causes  the  head  to  lose  energy,  expand,  and  erode.30 
The  estimate  given  here  shows  that  the  amount  of  beam  eroded  is  insignificant  relative  to  the 
pulse  duration. 


E.  Effective  Pinch  Current 

In  this  section  we  show  that  the  effective  pinch  current  is  essentially  equal  to  the  full  beam 
current  for  the  nominal  parameters  used  in  HEEB.  That  is,  the  plasma  formed  by  the  beam  in  air 
is  sufficiently  conducting  to  neutralize  the  beam  space  charge,  but  is  too  resistive  to  neutralize 
the  beam  current. 

A  beam  with  the  high  current  density  required  for  welding  and  materials  processing 
ionizes  and  heats  the  air  through  which  it  passes.  According  to  a  simple  ionization- 
recombination  model,  the  plasma  electron  density  ne  evolves  as 
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where  jb  =  V^a*,2  is  the  on-axis  beam  current  density,  de/dz  ■  3  keV/cm  is  the  collisional  energy 
loss  rate  discussed  in  Sec.  D-D,  W4  «  34  eV  is  the  average  energy  deposited  per  plasma  electron 
created,31  and  pT  ~  5xlO'8  cm3/s  is  the  dissociative  plasma  electron-ion  recombination 
coefficient.  After  a  characteristic  time 
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ne  on  axis  approaches  a  "saturated"  steady  state  value. 
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For  the  "welding"  beam  with  Ib  =  1  kA  and  aj,  =s  0.1  cm,  the  saturation  time  in  atmospheric  air  is 
ts  «  1  ns,  which  is  much  shorter  than  the  pulse  duration  rp  for  a  typical  case  with  rp  *  10  ns. 
The  plasma  density  ne,  and  conductivity  a  =  e^mv,  can  thus  be  treated  as  constant  over  most 
of  the  pulse;  here  v  is  the  plasma  electron  collision  frequency,  with  a  value  ~  2xl012  s'1  in 
atmospheric  air.  The  assumption  of  constant  a  is  used  later  to  simplify  the  analysis  of  hose 
instability.  In  broader  beams,  where  saturation  occurs  at  a  later  time,  it  may  be  necessary  to 
retain  the  time  dependence  of  a. 

The  plasma  neutralizes  the  beam  space  charge  in  a  time  l/4iro  which  is  negligibly  short 
(«  1  ns)  for  the  parameters  given  above.  Electrostatic  forces  can  therefore  be  safely  ignored. 
Current  neutralization,  however,  persists  only  for  a  time  t0  =  L/R  which  characterizes  magnetic 
field  diffusion  through  the  plasma.  Here  L  “  (l/c2)ln(l+a2/ab)  is  the  distributed  circuit 
inductance,  R  «  1/wa^o  is  the  distributed  plasma  resistance,  and  aw  >  ab  is  the  effective  outer 
edge32  of  the  plasma.  Since  R  steadily  decreases  as  beam  ionization  increases  the  plasma 
conductivity,  r0  is  an  increasing  function  of  the  retarded  time  r  =  t  -  z/Pc,  which  is  the  time 
measured  after  the  arrival  of  the  beam  head  at  location  z.  The  saturation  of  a  due  to 
recombination  eventually  limits  R,  and  this  restricts  r0  to 
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For  the  welding  beam  with  Ib  =1  kA,  ab  =  0.1  cm,  and  aw  <  50  ab,  the  decay  time  r0  is 
always  less  than  r.  The  plasma  return  current  induced  by  the  beam  then  decays  essentially 
instantaneously  so  that 


"ef  f 


(19; 


for  virtually  all  of  the  pulse.  All  but  the  early  beam  head  thus  feels  a  strong  magnetic  force  that 
keeps  the  beam  well  pinched  on  its  way  to  the  target.  For  broader  beams,  r0  is  larger  but  still 
smaller  than  the  pulse  duration  Tp  in  the  parameter  regime  of  interest.  In  that  case,  Ieff 
approaches  Ib  gradually  so  that  the  beam  pinch  force  increases  with  r.  This  gradual  increase 
causes  the  beam  head  to  flare,  so  that  ab  decreases  with  r.  Nordsieck  expansion  accentuates  the 
flaring  as  z  increases,  as  seen  in  the  simulations  presented  in  Sec.  IV. 
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in.  RESISTIVE  HOSE  INSTABILITY 


A.  General  Properties 

Electron  beams  propagating  through  dense  gas  are  subject  to  a  number  of  instabilities,  but 
the  only  significant  one  for  the  HEEB  applications  being  considered  here  is  the  resistive  hose 
instability.33'39  Resistive  hose  is  a  macro-instability  in  which  the  beam  performs  growing 
snakelike  oscillations  from  side  to  side,  in  response  to  magnetic  forces  from  the  plasma.  The 
instability  develops  as  plasma  eddy  currents  induced  by  the  beam  motion  lag  behind  the  beam  by 
a  magnetic  diffusion  time,  due  to  the  resistivity  of  the  plasma.  This  phase  lag  causes  the 
transverse  beam  wiggles  to  grow  in  amplitude.  Other  instabilities  are  known  to  occur  in  beams 
propagating  in  gas,  but  are  not  relevant  to  the  parameters  of  interest  here.  For  example,  the 
resistive  hollowing  instability  does  not  occur  in  the  absence  of  a  strong  plasma  return  current,40 
while  the  two-stream  instability  is  suppressed  by  collisions  between  the  plasma  electrons  and  the 
gas  molecules. 

Hose  instability  is  characterized  by  the  transverse  displacement  xb(z,T)  of  the  beam 
centroid  from  the  desired  propagation  axis.  When  xb  grows  to  large  amplitude  (xb  Jt  ab),  the 
instability  spreads  out  the  energy  density  delivered  to  the  target,  and  can  even  cause  the  beam  to 
miss  the  target.  To  avoid  these  effects,  the  beam  parameters  should  be  chosen  to  keep  hose  in  the 
linear  regime,  ^  <  a^  In  this  regime,  xb  is  proportional  to  the  initial  perturbation  Xj,  and  the  hose 
gain,  gh(z)  =  xmax(z)/Xj,  is  independent  of  the  magnitude  of  x(;  here  xmax(z)  is  the  maximum  of 
Xhfz.r)  within  the  pulse  at  a  given  propagation  distance  z.  Thus,  it  is  convenient  to  state  the  hose 
stability  requirement  as 

gh(Zt)  <  ^Xi  *  (20) 

Although  the  hose  gain  gh  is  independent  of  the  magnitude  of  x;  in  the  linear  regime,  it  is 
strongly  dependent35  on  the  temporal  variation  of  Xj  with  r.  To  be  conservative,  we  shall  allow 
for  excitation  of  the  most  rapidly  growing  hose  perturbations  and  compute  the  maximum  hose 
gain  possible.  Because  experience  shows  that  it  is  difficult  to  launch  a  beam  with  initial 
perturbations  less  than  a  few  percent  of  the  injected  radius  a0,  we  shall  restrict  gh  to  be  less  than 
30.  In  practice,  the  use  of  the  maximum  hose  gain  may  be  too  restrictive,  since  the  most 
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dangerous  hose  perturbations  usually  occur  at  frequencies  far  above  the  resonant  frequencies  of 
the  accelerator.  See  Sec.  III-E  for  further  discussion. 

In  the  next  three  subsections,  we  estimate  the  maximum  hose  gain  and  the  limitations  it 
imposes  on  the  choice  of  HEEB  parameters.  Particular  attention  is  given  to  three  stabilizing 
effects:  phase-mixing,  Nordsieck  expansion,  and  steering  by  a  metallic  pipe  surrounding  the 
beam.  The  latter  two  effects  are  shown  to  be  especially  important  for  applications  requiring 
narrow  beams. 


B.  Phase-Mixing  of  Electron  Trajectories 

Resistive  hose  develops  from  magnetic  interactions  between  the  beam  and  plasma  currents, 
each  of  which  is  displaced  from  the  axis  of  symmetry  and  from  each  other.  The  simplest  model 
of  the  instability,  which  we  shall  call  the  flexible-rod  model,  is  obtained  by  representing  the 
perturbations  as  displacements  xb(z,r)  of  the  beam  centroid,  with  no  internal  distortion  (and  thus 
no  phase  mixing),  and  characterizing  the  plasma  currents  by  the  location  of  the  net  current 
centroid  xn(z,r).  For  small  beam  displacements,  xb  «  ab,  a  linearized  equation  for  xb  can  be 
obtained  by  averaging  the  magnetic  force  over  the  radial  profile  of  beam  current.  In  the  absence 
of  wall  forces  (aw  »  ah),  this  yields35 
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where  ks  is  a  characteristic  hose  wavenumber  given  by 
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If  the  beam  and  plasma  have  identical  Bennett  profiles  given  by  Eq.  (4),  ks  =  2x(2l3)in/\fi, 
where  is  defined  in  Eq.  (5). 
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An  equation  for  the  net-current  centroid  xn  can  be  obtained  from  Ampere’s  law.  As 
derived  by  Lee35  (and  in  the  Appendix),  this  equation  is  given  in  the  limit  of  high  7  and  high  a 
by 


3x 

n 

Td  3  r 


(23) 


where  rd  is  a  magnetic  dipole  diffusion  time  given  by 

2  2 

Td  =  »Vo/2c 


(24) 


Here  a0  is  the  on-axis  conductivity.  Equations  (21),  (23),  and  (24)  are  appropriate  when  the 
plasma  return  current  is  negligible,  the  beam  and  plasma  have  identical  Bennett  profiles,  and 
xb  «  aj,. 

Equations  (21)  and  (23)  demonstrate  that  two  parameters  characterize  hose:  the  betatron 
wavelength,  «  k'1,  which  characterizes  the  response  of  the  beam;  and  the  dipole  diffusion 
time  Td,  which  characterizes  magnetic  relaxation  in  the  plasma.  In  general,  the  amount  of  hose 
growth  depends  on  the  propagation  distance  z,  scaled  to and  the  pulse  'ength  rp  scaled  to  rd. 
For  the  "welding"  beam  with  ab  =  0.1  cm,  Ib  =  1  kA  and  7  =  21,  we  find  X^,  *  12  cm  and 
rd  «  40  ps.  If  Zj  =  30  cm  and  rp  *  10  ns,  then  z/X^  «  2.5  while  Tp/rd  ^  250.  HEEB  welding 
thus  operates  in  a  long-pulse,  short-range  regime  where 


t  It  »  1, 

P  d 

(25a) 

2  /X_  *  few. 
t  p 

(25b) 

Materials  processing  applications  with  aj,  -  1  cm  also  typically  operate  in  this  regime. 

Most  of  the  previous  theoretical  work  on  hose  instability  has  been  in  the  opposite  short- 
pulse,  long-range  regime37  where 

r  It.  *  few,  (25c) 

P  d 

z  /Xa»  1.  (25d) 

t  p 
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Although  strategies  developed  to  control  hose  in  the  long-range  regime  are  usually  ineffective  in 
the  long-pulse  regime,  there  is  often  a  complementary  strategy  that  works.  The  general  rule  is 
that  variations  of  the  beam  parameters  in  t  are  important  when  (25c)  holds,  whereas  variations  in 
z  are  important  when  (25b)  holds.  For  example,  in  the  long-range  regime  (25d),  the  pulse 
duration  rp  must  be  limited  to  a  few  times  rd,  as  specified  in  (25c),  in  order  to  allow  the 
instability  to  convect35  out  the  beam  tail  before  reaching  intolerably  large  amplitude.  In  the 
long-pulse  regime  (25a)  of  interest  here,  it  is  the  propagation  distance  zt  that  must  be  limited,  to  a 
few  times  X^  as  specified  in  (25b).  In  a  similar  vein,  hose  growth  can  be  reduced  in  the  short- 
pulse  regime  (25c)  by  tailoring  the  beam  so  that  its  emittance,  and  therefore  its  radius  ab  and 
betatron  wavelength  X^,  decrease  significantly  (>  20%)  over  times  r  -  rd,  as  has  been  shown  by 
several  theoretical24,25,37  and  experimental26  studies.  In  the  long-pulse  regime,  however, 
tailoring  is  impractical  because  the  variations  needed  would  be  excessive  (as  would  be  the 
energy-deposition  radius  on  target).  Instead,  in  the  short-range  regime,  it  is  variations  in  X^  with 
z  (rather  than  t)  from  Nordsieck  expansion  that  helps  stabilize  hose.  A  last  distinction  between 
the  two  regimes  is  the  effect  of  Ib.  In  the  short-pulse  regime,  large  Ib  is  stabilizing  because  Td 
increases  with  Ib.  But  in  the  short-range  regime,  large  Ib  is  destabilizing  because  X^  decreases 
with  Ib.  In  this  paper,  we  concentrate  on  hose  growth  in  the  long-pulse,  short-range  regime. 

We  first  calculate  the  hose  growth  rate  using  Eqs.  (21)-(23),  which  omit  phase-mixing  of 
the  electron  trajectories.  Since  the  beam  and  plasma  equilibria  are  assumed  to  be  independent  of 
z  and  r,  we  consider  solutions  of  the  form  x  **  e'lkz‘l£JT.  We  take  u>  to  be  real  and  calculate  the 
corresponding  complex  value  of  k,  because  our  primary  interest  is  to  calculate  hose  growth  in  z. 
Equations  (21)  and  (23)  yield  the  flexible-rod  dispersion  relation, 


2  1WTd 

(k/k  )=  -  ^ - 

s  l-i ur 


(26) 


Writing  k  in  terms  of  real  and  imaginary  components,  k  =  kr+ik(,  we  find  that  the  maximum 
value  of  kj  is  given  by 

kmaX  =  k  /f8  .  (27) 

1  s 
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This  defines  the  peak  growth  rate  in  z,  in  the  absence  of  phase-mixing  of  the  electron  trajectories. 

Phase-mixing  is  a  particle  effect  that  occurs  because  the  beam  pinch  force  is  anharmonic 
(i.e.,  nonlinear  in  r),  which  makes  the  betatron  wavelength  of  each  beam  electron  depend  on  the 
amplitude  of  its  orbit.  The  spread  of  betatron  wavelengths  for  the  electrons  in  a  given  beam  slice 
destroys  the  coherence  of  hose  oscillations  in  that  slice,  thereby  damping  those  oscillations. 
Stated  differently,  phase-mixing  gradually  converts  the  ordered  energy  of  hose  motion  into  beam 
thermal  energy  and  emittance  (with  the  increase  in  emittance  being  small  when  xb  «  ah). 

To  incorporate  phase-mixing  in  an  analytically  tractable  way,  we  use  the  spread-mass  model 
developed  by  Lee.35  In  this  model,  the  beam  is  artificially  represented  as  a  superposition  of 
"flexible-rod"  components  corresponding  to  different  electron  masses,  and  thus  different  betatron 
wavelengths.  The  distribution  of  electron  masses  is  chosen  to  preserve  the  average  oscillation 
frequency  ks.  This  model  has  been  shown  to  give  hose  growth  rates  in  good  agreement  (~  10 
percent)  with  those  from  more  accurate  but  involved  treatments.38  An  abbreviated  description  of 
the  model  is  given  in  the  Appendix.  Here  we  simply  quote  Lee’s  result  [Eq.  ( 143)  of  Ref.  35]: 

kmaX  =  0.183  kfl  =  0.317  k  ,  (28) 

i  Po  s 

where  kpo  is  the  maximum  (on-axis)  betatron  oscillation  frequency.  This  peak  rate  is  roughly  ten 
percent  smaller  than  prediction  (27),  a  difference  comparable  to  the  accuracy  of  the  spread-mass 
model  itself.  Thus,  phase-mixing  reduces,  but  only  slightly,  the  hose  growth  rate  in  z. 

Phase-mixing  has  a  weak  effect  on  hose  growth  in  z  because  it  introduces  damping  at  a  rate 
that  is  small  relative  to  ks.  The  damping  has  a  pronounced  effect,  however,  on  growth  in  r  :  by 
spreading  the  resonance  at  ks,  it  makes  the  instability  convective  rather  than  absolute.35  This 
effect  is  critical  in  the  long-range  regime,  but  not  in  the  long-pulse,  short-range  regime  of  interest 
here. 
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C.  Nordsieck  Expansion 

The  preceding  analysis  shows  that  hose  growth  increases  exponentially  with  the  number  of 
betatron  wavelengths  X^  propagated.  As  a  result,  any  process  that  increases  \e  suppresses  hose. 
One  such  process  is  Nordsieck  expansion  which  causes  the  beam  radius  ab,  and  therefore  X^ 
(which  is  proportional  to  a^),  to  increase  with  propagation  distance  z. 

A  simple  estimate  of  the  effects  of  Nordsieck  expansion  on  hose  can  be  easily  made  in  the 
"welding"  regime.  In  that  regime,  7  and  Ieff  are  nearly  constant  with  both  z  and  r,  and  thus  the 
Nordsieck  length,  LN  =  LN(7leff),  can  be  taken  as  constant.  From  Eqs.  (5)  and  (9),  the  betatron 
wavelength  can  be  taken  as  independent  of  r  but  growing  in  z  as 

X  (z)  =  X  -  (0)  exp  (z/L  )  .  (29) 

P  P  N 

If  hose  grows  nearly  exponentially  with  the  number  of  X^  propagated,  we  expect  the  hose  gain 
ghfz,)  to  be  given  by 

f*‘  dz 

ln[gh(ztH  -  «h  I  J~u7 

0  p 

=  xrrfr  l1  -  exp(-w]  •  ,30) 

P 

where  «h  is  a  hose  growth  coefficient  of  order  unity,  whose  precise  form  will  be  discussed  in  the 
next  two  subsections  and  specified  in  Eq.  (39). 

Equation  (30)  indicates  that  Nordsieck  expansion  begins  to  reduce  hose  significantly  at 
distances  zt  *  LN/5,  and  that  it  limits  the  hose  gain  to  gh  s  exp[ahLN/X^(0)].  Although  a  more 
formal  analysis37  shows  that  Eq.  (30)  is  strictly  valid  only  when  LN  >  zt  and  LN  »  \p(0),  the 
restriction  is  usually  unimportant  since  hose  growth  is  modest  outside  that  regime. 

D.  Wall  Stabilization  by  a  Metallic  Pipe 

A  metallic  pipe  or  "straw"  enclosing  the  beam  further  suppresses  hose  as  return  currents 
induced  in  the  pipe  walls  magnetically  steer  the  beam.  The  return  current  in  the  walls  of  the 
straw  can  be  represented  in  the  linearized  regime  as  an  image  line  current  of  magnitude  -I„ 
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located  at  a^/xn  from  the  center  of  the  straw.  Here  the  net  current  In  equals  the  sum  of  beam  and 
plasma  return  currents,  and  aw  >  aj, »  xn  is  the  radius  of  the  straw.  In  the  present  application, 
we  have  assumed  that  In  «  Ib. 

The  return  current  in  the  straw  adds  another  force  component  to  Eq.  (21): 


*b 


3z' 


=  k  (x  — x,  )  - 
s  no 


'b  v2 
-r  k  x 

2  s  n 


w 


(31) 


where  rb  is  a  characteristic  beam  radius  defined  by 


(32) 


In  general,  rb  is  comparable  to  but  not  equal  to  the  half-current  radius  ah.  If  the  beam  and  plasma 
have  identical  Bennett  profiles  of  radius  a,,,  then  rb  =  'TTab. 

A  second  effect  of  the  straw  is  to  accelerate  the  diffusion  of  the  dipole  magnetic  fields.  As 
shown  in  the  Appendix,  the  magnetic  dipole  decay  time  is  reduced  to 

Td‘  [1'(ab/aw)2ln(1*aw/ab,](’rab<,o/2c2)  ‘  <33) 

in  place  of  Eq.  (24).  However,  since  rd  was  already  much  smaller  than  the  pulse  duration  rp,  the 
reduction  has  no  effect  on  the  hose  gain  gh. 

For  simplicity,  in  the  present  discussion  we  ignore  phase-mixing  (included  in  the  Appendix) 
and  Nordsieck  expansion.  Again  considering  solutions  of  the  form  x  ~  e',kz‘,WT  with  real  u  and 
complex  k,  Eqs.  (21)  and  (31)  give  a  flexible-rod  dispersion  relation. 


2  2 

o  1-r.  /a* 

(k/k  )2  =  1  -  .b  W 
S  1-1CJT  , 


(34) 


in  place  of  Eq.  (26).  The  peak  growth  rate  is  then  computed  to  be 
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(35) 
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2  2 
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2  2 
1+r^/a 
b  w 


in  place  of  Eq.  (27).  Equation  (35)  indicates  that  wall  forces  reduce  the  hose  growth  rate  for 
Bennett  profiles  by  a  factor 


,  max  . 
k.  (a  ) 
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(36) 


This  reduction  appears  to  be  consistent  with  that  given  by  Uhm  and  Lampe39  for  flat-topped 
beams. 

In  Fig.  4  the  dashed  curve  shows  the  normalized  peak  growth  rate  k^/k^  from  Eqs.  (27) 
and  (36)  for  a  Bennett  beam  as  a  function  of  a^/a^,,  where  kjD  =  TTk,.  The  solid  curve  shows  the 
corresponding  result  as  calculated  in  the  Appendix,  with  phase  mixing  included.  The  two  results 
are  in  good  qualitative  agreement,  indicating  that  the  effect  of  phase  mixing  is  small.  In  a  typical 
implementation  with  a^,  =  38^  Fig.  4  indicates  that  the  presence  of  the  straw  cuts  the  peak  growth 
rate  by  nearly  a  factor  of  two.  This  allows  the  beam  to  propagate  nearly  twice  as  far  in  z  for  a 
given  hose  gain  gh,  an  improvement  that  far  exceeds  the  stabilizing  effect  of  either  phase-mixing 
or  Nordsieck  expansion. 

In  addition  to  improving  stability,  the  straw  could  be  used  for  a  variety  of  other  purposes, 
including  the  following:  (i)  Scattering  can  be  reduced  by  flowing  helium  through  the  straw,  or 
by  differentially  pumping  to  reduce  the  air  pressure  by  a  factor  of  five  to  ten.  As  we  shall  see, 
this  can  substantially  reduce  the  required  beam  electron  energy.  Flowing  helium  through  a  straw 
should  be  simpler  and  more  economical  than  blowing  a  helium  jet2  across  open  air.  (ii)  The 
straw  can  be  used  to  shield  the  beam  from  stray  magnetic  fields4  as  discussed  in  Sec.  V-B.  (iii) 
The  straw  can  be  used  for  precision  corrections  to  beam  pointing.  A  straight  straw  can  accept  a 
misaligned  beam  and  magnetically  steer  it  through  a  (typically  small)  angle 


dx/dz  ^ 


<  W  <2VV 


1/2 


(37) 
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Alternatively,  a  curved  straw  can  direct  the  beam  around  a  radius  of  curvature  which  satisfies 

R  *  (a2/a,  )  (I _ /2I.  )  .  (38) 

w  b  A  b 

Equations  (37)  and  (38)  were  derived  by  requiring  that  the  beam  remain  within  one  beam  radius 
of  the  straw  axis. 

E.  Estimating  the  Hose  Gain 

The  results  of  the  three  previous  sections  on  hose  may  be  combined  by  setting  the  hose 
growth  coefficient  ah  in  Eq.  (30)  to 

(39) 

b  w 

where  1.63  =  2ir(0.317)(2/3)1/2.  This  value  was  obtained  using  the  spread-mass  result  (28)  for 
Id-1"1  in  terms  of  ks,  and  using  Eq.  (36)  for  the  reduction  in  hose  growth  produced  by  a  straw  of 
radius  aw,  assuming  Bennett  profiles.  Equations  (30)  and  (39)  allow  us  to  estimate  the  hose  gain 
gh  as  a  function  of  the  propagation  distance  zt,  the  initial  betatron  wavelength  X^(0)  given  by 
Eq.  (5),  and  the  Nordsieck  length  LN  given  by  Eq.  (10). 

For  the  canonical  welding  parameters  of  z,  =  30  cm,  ^  =  1  kA,  7  =  21,  and  a0  =  0. 1  cm,  the 
predicted  hose  gain  is  gh  <  35  in  open  air,  and  it  drops  to  gh  <  8  if  a  straw  of  radius 
aw  =  3ab  =  0.3  cm  is  used.  In  the  former  case,  the  transverse  beam  displacements  Xj  at  the 
accelerator  would  have  to  be  less  than  a  few  percent  of  a0  to  satisfy  criterion  (20)  at  the  target;  in 
the  latter  case,  Xj  could  be  as  large  as  10  percent  of  a0,  which  is  well  within  present  accelerator 
capabilities.  Moreover,  as  mentioned  in  Sec.  III-A,  the  restrictions  on  x(  are  conservative, 
because  they  apply  only  if  the  most  rapidly  growing  hose  frequencies  in  r  are  excited.  For  the 
canonical  welding  beam,  the  most  dangerous  hose  frequencies  (wrd  -  1)  occur  at  a  few  GHz.  If 
these  frequencies  are  not  excited  or  are  excited  only  at  very  low  levels,  the  restrictions  on  x( 
could  be  relaxed. 
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IV.  PROPAGATION  EXAMPLES 


A.  HEEB  Parameter  Choice 

As  discussed  earlier,  the  two  main  propagation  issues  affecting  HEEB  operation  are  beam 
expansion  and  the  resistive  hose  instability.  Expansion  must  be  limited  to  prevent  the  beam  from 
becoming  too  broad  on  target,  and  hose  growth  must  be  limited  to  prevent  the  beam  from 
wandering  off  target.  We  have  expressed  the  former  limit  in  terms  of  an  allowed  expansion 
factor  ga  and  the  latter  limit  in  terms  of  an  allowed  hose  gain  gh.  As  Eqs.  (5),  (12)  and  (30) 
indicate,  self-pinched  expansion  restricts  7leff  while  the  hose  requirement  restricts  Tr/Ieff.  In 
addition,  we  limit  the  beam  energy  to  y  a  21,  in  order  to  avoid  nuclear  activation  of  the  target. 
Together,  these  three  limits  circumscribe  a  range  of  values  allowed  for  7  and  Ieff. 

The  range  of  allowed  parameters  is  illustrated  in  Fig.  5  which  shows  three  curves 
representing  the  limits  imposed  by  beam  expansion  (solid  curve),  resistive  hose  (dashed  curve), 
and  avoidance  of  nuclear  activation  (dot-dashed  curve).  In  this  parameter  regime,  Ieff « Ib,  so  that 
the  results  are  plotted  in  terms  of  Ib  versus  beam  energy.  In  all  of  the  cases  shown  in  Fig.  5,  it  is 
assumed  that  the  beam  is  injected  in  equilibrium  with  a  radius  of  aG  =  0.1  cm,  that  the  expansion 
on  target  is  restricted  to  ga  S  FT,  and  that  propagation  is  in  full-density  air.  In  Fig.  5a,  the 
propagation  distance  is  Zj  =  10  cm,  the  allowed  hose  growth  is  gh  £  30,  and  there  is  no  straw.  In 
that  case  the  (shaded)  area  circumscribed  by  the  three  curves  is  large,  indicating  considerable 
flexibility  in  the  choice  of  beam  current  and  energy.  In  Fig.  5b,  the  propagation  distance  is 
increased  to  z,  =  20  cm,  which  reduces  the  allowed  parameter  range.  At  z,  =  30  cm  (Fig.  5c),  the 
beam  parameters  are  tightly  constrained  to  Ib  •*  1  kA  and  e  «  10  MeV. 

A  wide  parameter  range  for  propagation  to  zt  =  30  cm  can  be  recovered  by  propagating  the 
beam  inside  a  conducting  straw  of  radius  aw  =  3  mm.  As  shown  in  Fig.  6a,  the  beam  energy  can 
then  be  as  low  as  e  *  6  MeV,  with  the  current  increased  to  Ib  *  1 .6  ka.  With  the  use  of  a  straw  to 
guide  the  beam,  a  reasonable  choice  of  parameters  is  available  even  if  the  criterion  for  acceptable 
hose  growth  is  reduced  to  a  conservative  gh  £  10.  The  allowed  parameter  range  is  substantially 
larger  still  if  the  straw  is  Filled  with  helium,  or  with  air  at  reduced  density.  The  case  of  He  at 
atmospheric  density  (equivalent  to  air  at  40torr)  is  shown  in  Fig.  6b,  which  indicates  that  the 
beam  energy  can  be  as  low  as  1  MeV.  In  addition,  in  all  cases  the  required  propagation 
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parameters  can  be  relaxed  if  the  beam  radius  is  allowed  to  be  several  mm  or  more  on  target.  (In 
fact,  the  on-target  radius  may  have  to  be  larger  than  1  mm,  in  order  to  avoid  shocking  the 
material,  for  typical  beam  pulse  lengths  of  several  tens  of  ns;  moreover,  radii  ~  2  mm  do  not 
appreciably  compromise  weld  strength,  but  do  reduce  the  need  for  tight  joint  tolerances.) 
Overall,  we  conclude  that  narrow  beams  of  high  energy  density  can  be  readily  delivered  to 
targets  located  30  cm  or  more  from  the  accelerator,  without  need  of  vacuum  pumping. 

B.  Particle-Simulation  Results 

The  analytical  estimates  given  for  beam  expansion  and  hose  instability  in  the  previous 
sections  make  a  number  of  simplifying  assumptions  which  are  marginal  in  some  cases.  To  relax 
these  assumptions  and  obtain  a  more  detailed  treatment,  we  use  the  three-dimensional,  nonlinear 
particle  simulation  code  SARLAC,  which  incorporates  gas  scattering,  ohmic  and  collisional 
energy  loss,  and  a  set  of  rate  equations  for  computing  the  plasma  conductivity.41  The  beam  is 
taken  to  be  ultrarelativistic  (y  »  1)  and  paraxial  (v±  «  vz),  so  the  retarded  time  t  =  t  -  z/c  is  a 
constant  of  the  motion  for  the  beam  particles.  The  beam  propagates  in  a  gas-filled  conducting 
straw  with  a  radius  aw  that  can  be  much  larger  than  the  beam  radius  ab.  The  independent 
variables  are  r,  d,  r,  and  z,  where  $  is  the  azimuthal  angle  about  the  axis  of  the  straw.  With 
these  assumptions.  Maxwell’s  equations  and  the  conductivity  rate  equations  contain  no 
derivatives  in  z,  whereas  the  equations  of  motion  for  the  beam  particles  contain  no  derivatives  in 
t.  The  code  allows  for  variations  in  the  beam  radius,  current,  and  energy  within  the  pulse. 

The  beam  was  injected  with  current  Ib  rising  rapidly  to  1  kA,  energy  7  =  21,  radius 
aD  =  0.1  cm  (but  slightly  larger  at  early  t),  and  pulse  length  shortened  to  Tp  =  2  ns  to  conserve 
computer  time.  The  beam  was  loaded  with  a  Bennett  radial  distribution,  and  a  transverse 
temperature  matched  to  pinch  current  Ieff  =  0.85  kA.  At  z  =  0,  initial  sinusoidal  hose 
displacements  Xj(r)  and  y{(r)  were  applied  in  both  the  x  and  y  directions.  The  amplitude  of  these 
initial  perturbations  was  10'3  cm,  and  the  frequency  of  x{  was  set  to  cj/2k  =  1  GHz,  while  the 
frequency  of  y;  was  set  to  5  GHz.  (In  the  linear  regime,  the  perturbations  of  the  beam  in  the  x-z 
plane  and  in  the  y-z  plane  evolve  independently.)  The  use  of  a  shortened  rp  and  only  two 
discrete  hose  perturbation  frequencies  leads  to  a  slight  underestimate  of  the  maximum  possible 
hose  growth. 
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In  the  first  set  of  simulations,  the  radius  of  the  straw  was  set  to  aw  =  5  cm,  which  is  large 
enough  so  that  wall  effects  should  be  negligible.  The  beam  half-current  radius  a^T)  is  plotted  at 
two  z  locations  in  Fig.  7a.  As  is  evident,  the  beam  slowly  expands  and  its  head  gradually  erodes 
as  it  propagates.  Hose  growth  for  the  slice  at  r  =  0.75  ns  is  shown  by  the  solid  curve  in  Fig.  7b. 
Maximum  hose  growth  occurred  at  somewhat  later  r,  with  a  peak  growth  of  gh  =  25  at 
z  =  30  cm.  This  agrees  well  with  the  analytic  prediction  of  gh  <  35  given  earlier  (for  the  most 
unstable  perturbation  frequency  and  unlimited  pulse  length). 

In  the  second  set  of  simulations,  the  same  beam  [with  its  density  nb(r)  cut  off  at  the  straw 
radius  aw]  was  propagated  inside  a  conducting  straw  of  radius  aw  =  3  mm.  As  shown  by  the 
dashed  curve  in  Fig.  7b,  the  straw  reduced  the  peak  hose  growth  to  gh  =  5  at  z  =  30  cm.  Again, 
this  agrees  well  with  the  analytic  prediction  of  gh  <  8.  The  loss  of  beam  current  due  to  scrape-off 
on  the  straw  is  shown  in  Fig.  7c  to  be  only  five  percent. 

To  test  for  wall  steering,  the  preceding  simulation  was  repeated  except  that  the  beam  was 
injected  misaimed  at  an  angle  dx/dz  =  10  mrad.  As  shown  in  Fig.  7d,  the  beam  slice  at 
t  =  0.75  ns  oscillated  about  a  peak  displacement  of  0.35  mm  with  a  wavelength  X  «  28  cm.  This 
displacement  and  wavelength  agree  well  with  analytic  prediction.  The  figure  also  shows  that 
hose  growth  is  unaffected  by  the  misaiming. 
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V.  OTHER  ISSUES 


A.  Gas  Heating 

In  gas,  electron  beams  deposit  energy  directly  into  excited  states  and  the  plasma  electron 
population,  and  this  energy  is  shared  later  with  the  translational,  vibrational,  and  other  degrees  of 
freedom  of  the  gas.  The  heating  and  equilibration  of  the  gas  produces  a  reduced-density  channel, 
with  an  eventually  elevated  plasma  conductivity  a.  For  a  single  isolated  HEEB  pulse  (which 
may  be  suitable  for  some  applications),  these  effects  have  no  significant  impact  on  propagation, 
because  the  pulse  length  is  too  short  for  appreciable  gas  expansion  to  occur,  and  because 
increases  in  o  in  the  beam  tail  have  virtually  no  effect  on  pinching  or  overall  hose  growth. 

However,  in  a  repetitive  mode,  late  beam  pulses  encounter  the  hot  air  channel  formed  by 
their  predecessors,  and  the  effect  of  this  channel  on  beam  propagation  depends  sensitively  on  the 
channel  temperature.  At  channel  temperatures  below  3000  K  (and  pulse  repetition  rates  below 
10  kHz),  the  plasma  recombines  between  pulses,  and  the  heated  channel  expands.  The  beam 
pulses  then  enter  a  rarefied  but  non-conducting  channel  that  retards  Nordsieck  expansion  of  the 
beam.  In  addition,  the  radial  density  gradient  in  the  channel  alters  the  beam-generated 
conductivity,  and  this  leads  to  a  weak  channel-guiding  force  on  the  beam.42,43 

At  temperatures  above  -  3500  K,  the  effect  of  the  channel  changes  abruptly.  At  these 
temperatures,  the  channel  remains  ionized  and  highly  conducting  between  pulses.  A  fresh  beam 
pulse  entering  the  preionized  channel  induces  a  plasma  return  current  that  fully  neutralizes  the 
beam  current,  and  this  current  magnetically  ejects  the  beam  from  the  channel.  The  transition 
from  channel  guiding  to  channel  ejection  has  been  observed43  to  occur  at  a  temperature 
-  3500  K.  Very  hot  channels  thus  deflect  subsequent  beam  pulses  away  from  the  target.  At  the 
high  current  densities  required  for  HEEB  welding,  channel  temperatures  can  reach  3500  K  after  a 
few  pulses. 

One  way  to  avoid  channel  overheating  is  to  blow  the  channel  away.  An  unlimited  number 
of  beam  pulses  can  then  propagate,  since  each  pulse  encounters  essentially  virgin  gas.  An  air 
flow  speed  of  a  few  m/s  across  the  channel  is  sufficient  for  pulse  rates  less  than  a  few  kHz.  In 
open  air,  natural  convection  may  provide  sufficient  cooling.  If  the  beam  is  propagated  in  a  straw, 
it  may  be  necessary  to  blow  gas  through  the  straw,  or  out  through  perforations. 
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B.  Beam  Deflection 

With  the  low-energy  beams  presently  used  for  welding,  beam  deflection  by  stray  magnetic 
fields  is  a  significant  problem.4,6  8  The  high-energy  beams  discussed  in  this  paper  are  so  stiff 
that  they  are  virtually  impervious  to  naturally  occurring  magnetic  fields  localized  on  the 
workpiece  surface.  However,  magnetic  fields  induced  by  the  beam,  and  external  fields  that  act 
coherently  over  the  entire  beam  path,  may  still  be  of  significance.  A  convenient  measure  of  these 
effects  is  the  value  of  the  perpendicular  magnetic  field  B0  which,  acting  over  the  entire 
propagation  length,  is  sufficient  to  deflect  the  beam  by  an  amount  equal  to  its  radius  ab: 
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For  our  canonical  HEEB  parameters,  B0  =  8  G.  The  earth’s  magnetic  field  (0.3  G)  is  therefore 
unimportant,  and  deflection  from  power  lines  is  usually  unimportant  provided  the  cables  lie  a 
meter  or  more  away  from  the  beam. 

A  more  significant  concern  is  the  image  currents  and  charges  induced  by  the  beam  in  nearby 
conductors.  Consider  a  beam  propagating  in  open  air  parallel  to  a  ground  plane.  The  net  current 
(the  sum  of  Ib  and  the  induced  plasma  current)  induces  an  equal  but  opposite  image  current  in  the 
ground  plane,  and  this  current  produces  a  magnetic  field  of  strength  Bg  =  In/xc  on  the  beam 
located  a  distance  x  away.  This  magnetic  field  would  tend  to  deflect  the  beam  away  from  the 
ground  plane.  The  field  Bg  exceeds  B0  when  x  s  (zt2/2ab)(In/IA).  For  our  canonical  welding 
beam  parameters,  this  occurs  when  x  S  12.6  cm.  Thus  this  simple  calculation  would  appear  to 
indicate  that  magnetic  deflection  could  be  a  significant  problem. 

Fortunately,  there  are  a  number  of  effects  that  can  greatly  reduce  the  strength  of  the 
interaction  between  a  beam  and  a  nearby  ground  plane.  A  brief  catalog  of  these  effects  would 
include  the  following:  (i)  Image  charges  are  also  induced  in  the  conductor,  and  the  resulting 
electrostatic  forces  are  opposite  in  sign  to  the  magnetic  image  forces.  If  the  beam  is  highly 
relativistic,  these  forces  are  very  nearly  equal  and  opposite  to  the  magnetic  forces,  (ii)  However, 
the  plasma  channel  surrounding  the  beam  can  shield,  to  differing  degrees,  both  the  magnetic  and 
electrostatic  image  forces.  Usually,  the  electrostatic  forces  are  shielded  much  more  effectively 
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than  the  magnetic  forces,  (iii)  Under  some  circumstances,  there  are  significant  time  delays  for 
the  transmission  of  fields  from  the  beam  to  the  conductor  and  back,  (iv)  The  deflecting  force  on 
the  beam  rises  from  zero  with  r,  due  to  variations  in  Ib(r),  o(t),  etc.  Thus  the  deflecting  forces 
tend  to  send  different  parts  of  the  beam  off  at  different  trajectories,  (v)  However,  there  are  also 
strong  coupling  forces  (mediated  by  the  plasma  channel)  that  hold  the  beam  together,  and  force 
the  body  of  the  beam  to  follow  the  trajectory  established  by  the  head.  In  a  recent  paper,44  we 
sorted  through  these  issues  in  considerable  detail,  and  concluded  that  ground  plane  deflection  is 
essentially  absent  when  two  conditions  are  met:  the  beam  must  be  highly  relativistic,  7  >  5;  and 
it  must  have  sufficient  current  to  ionize  the  gas  rapidly,  in  the  sense  that 

s.I  >  a  /x  ,  (41) 

ib  b 

where  S;  =  (e/2mvWjC2)(de/dz)  is  a  beam  ionization  coefficient.  Under  these  conditions,  the 
beam  body  follows  the  trajectory  determined  by  the  head.  At  the  head  of  the  beam,  there  is  not 
enough  plasma  to  neutralize  the  beam  fields,  and  thus  the  head  is  not  deflected  significantly 
because  the  electric  and  magnetic  image  forces  cancel  (to  order  7'2).  In  air,  Sj  «  0.04  kA'1,  which 
leads  to  the  prediction  that  ground-plane  deflection  becomes  unimportant  at  distances  x  >  3  cm 
for  our  canonical  welding  beam.  If  it  is  necessary  to  propagate  the  beam  closer  than  3  cm  to  a 
ground  plane,  the  beam  should  be  propagated  inside  a  conducting  straw,  to  shield  it  from  outside 
induced  fields. 

An  additional  deflection  mechanism  arises  in  welding  from  asymmetrical  currents  induced 
in  the  melt  pool.8  If  these  currents  were  confined  to  one  side  of  the  joint,  the  beam  could  be 
deflected  by  its  radius  in  a  penetration  distance  as  short  as  X^/4,  which  is  3  cm  for  the  nominal 
parameters  given  earlier.  In  practice,  however,  the  distance  is  much  longer  because  the  currents 
are  not  unilaterally  distributed.  Consequently,  beam  deflection  within  the  melt  pool  should  not 
be  a  problem,  even  for  welds  as  deep  as  10  cm. 

C.  Matching 

In  the  preceding  calculations,  we  assumed  that  the  beam  enters  in  pinched  equilibrium  with 
Tx  =  Tb,  a  Bennett  radial  profile,  and  a  radius  suitable  for  the  application.  If  the  beam  is  injected 
out  of  equilibrium,  the  beam  radius  and  profile  will  oscillate.  However,  as  described  in  the  hose 
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analysis,  these  oscillations  will  damp  as  the  beam  electrons  get  out  of  phase  with  one  another  in 
the  anharmonic  pinch  potential.  Phase-mixing  thus  causes  the  beam  to  approach  a  quiescent 
steady  state  after  roughly  one  betatron  wavelength  of  propagation.45  Phase-mix  damping 
increases,  however,  the  beam  emittance.24 

A  beam  injected  with  a  matched  temperature  T±  =  TB  but  a  non-equilibrium  radial  profile 
will  experience  only  modest  emittance  growth,45  provided  the  plasma  return  current  is  much 
smaller  than  Ib,  as  is  true  for  the  parameters  considered  here,  and  provided  the  injected  profile 
falls  off  smoothly  with  r.  Thus  it  is  not  particularly  important  for  our  applications  to  launch  a 
beam  whose  radial  profile  is  well-matched. 

On  the  other  hand,  to  prevent  the  beam  from  over-expanding  at  the  target,  the  temperature  at 
injection  should  kept  at  or  below  TB.  If  the  beam  is  injected  cold,  T±  <  TB>  it  will  collapse  under 
the  action  of  the  pinch  force,  but  the  transverse  energy  acquired  during  collapse  will  push  the 
electrons  back  out,  so  that  the  equilibrated  radius24  is  smaller  than  the  injected  radius  only  by  a 
factor  ~  TT.  That  is,  the  beam  acquires  enough  emittance  during  equilibration  to  prevent 
significant,  permanent  collapse.  A  related  point  is  that  the  injected  radius  a^,  should  be  no  larger 
than  that  desired  at  the  target,  since  emittance  growth  prevents  the  beam  from  contracting  much 
below  a,,.  At  the  same  time,  a^,  should  not  be  much  smaller  than  needed,  since  this  increases  hose 
growth.  In  general,  therefore,  the  beam  should  be  injected  with  a„  slightly  less  than  needed  at  the 
target,  and  with  T±  £  TB. 
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VL  CONCLUSION 

The  advent  of  high-power,  high-energy,  repetitive  electron  beams  offers  new  possibilities 
for  electron  beam  welding  and  materials  processing.  In  this  report,  we  have  focused  on  the 
potential  of  these  beams  to  propagate  tightly  pinched  (with  radius  -  1  mm)  through  a  few  tens  of 
cm  in  ambient  air,  without  need  of  external  guide  fields.  This  capability  would  eliminate  the 
need  for  vacuum  pumping  at  the  target  and  for  magnetic  focusing  at  the  accelerator. 
Furthermore,  the  use  of  high-energy  electrons  offers  the  potential  for  improved  deep  welding, 
due  to  the  long  range  of  the  electrons  in  the  target  material  and  the  constancy  of  the  beam  radius 
in  the  "keyhole"  made  within  the  target. 

After  reviewing  the  physics  of  self-pinched  propagation  in  gas,  our  major  conclusion  is  that 
propagation  is  an  important  but  not  overriding  issue  for  most  applications.  That  is,  the 
limitations  imposed  by  propagation  are  significant  but  do  not  severely  restrict  the  utilization  of 
the  technology.  The  three  most  important  limitations  are:  gas  scattering,  which  limits  the 
minimum  beam  power;  hose  instability,  which  limits  the  maximum  beam  current;  and  the  need  to 
avoid  nuclear  activation,  which  limits  the  maximum  beam  energy.  The  restrictions  on  beam 
current  and  energy  become  stressing  if  the  beam  is  required  to  remain  small  (aj,  ~  0. 1  cm)  while 
propagating  over  a  significant  standoff  to  the  target.  The  nominal  parameters  needed  to 
propagate  through  30  cm  of  full  density  air  are  1  kA  and  10  MeV;  these  values  are  sensitive  to 
the  desired  beam  range  and  can  be  reduced  if  the  standoff  distance  is  reduced.  The  margin  of 
safety  can  be  increased  by  propagating  the  beam  inside  a  conducting  "straw"  which  reduces  hose 
growth,  eliminates  concerns  regarding  beam  deflection  from  stray  or  beam-induced  fields,  and 
can  reduce  gas  scattering  if  the  straw  is  differentially  pumped  or  filled  with  a  gas  like  helium.  In 
the  regime  of  interest,  the  beam  pulse  length  is  not  limited  by  propagation  considerations. 

In  short,  the  advantages  of  high-power  electron  beams  for  welding  and  materials  processing 
can  apparently  be  realized,  provided  pulsed  high-power  beam  sources  become  available  with  cost 
and  reliability  comparable  to  that  of  the  low-power  dc  generators  being  used  commercially  today. 
These  advantages  include:  the  ability  to  propagate  to  the  target  in  ambient  air  without  need  of 
vacuum  pumping  or  magnetic  focusing;  the  ability  to  deposit  energy  deep  within  the  target;  the 
ability  to  heat  the  target  rapidly  to  minimize  thermal  losses  and  to  excite  shocks  if  desired;  and 
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the  ability  to  fine-tune  the  time  history  of  the  deposition.  Finally,  it  is  quite  feasible  to  test  the 
propagation  predictions  made  here  using  scaled  experiments  which  can  be  performed  with 
presently  available,  high-power  pulsed  beam  sources. 
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APPENDIX:  RESISTIVE  HOSE  ISSUES 


In  th's  appendix,  we  refine  the  analysis  given  for  the  resistive  hose  instability  in  the 
presence  of  a  conducting  straw.  In  particular,  we  show  that  the  straw  reduces  the  dipole  decay 
time  rd  according  to  Eq.  (33),  and  we  include  the  effects  of  orbit  phase-mixing. 

To  compute  the  dipole  decay  time,  we  drop  electrostatic  effects  and  axial  derivatives  from 
Ampere’s  law,  leaving35 


A 

z 


3A 

a  z 
c  dr 


(Al) 


where  Az  is  the  axial  vector  potential.  The  neglect  of  electrostatics  is  justified  provided  the 
plasma  conductivity  is  high,  4»a0ab/c  »  1,  while  the  neglect  of  axial  derivatives  is  justified 
provided  the  beam  is  energetic,  7  »  1,  and  paraxial,  »  a,,.  To  simplify  the  analysis  further, 
we  assume  that  the  plasma  conductivity  a  follows  the  beam  and  is  distributed  with  the  same 
radial  profile.  The  straw  imposes  the  boundary  condition  that  =  0  at  r  =  aw. 

For  small  displacements  xb  from  the  axis  of  the  straw,  the  beam  current  density  can  be 
expanded  in  the  (r,0)  plane  as 


j  (r,0fr)  =  jb<r,0,T)  -  x^t)  cos 9  jb(r,0,r)  + 


(A2 ) 


The  plasma  conductivity  a  can  be  expanded  similarly,  while  the  potential  Az  can  be 
expanded  as 

3A 

A  (r ,0,t)  =  A  (r ,9,t)  +  A  (r,T)  -  x  (r)  cos0  -r— —  +  ...  (A3) 

z  w  o  n  or 


Here  xn  is  the  displacement  of  the  net-current  centroid,  and  Aw  is  the  contribution  to  the  vector 
potential  due  to  the  image  current  -I,,  located  at  ^2/xn  from  the  axis  of  the  straw.  For  xn  «  aw 
and  r  <  aw, 

21 

A  (r ,B,t)  =  -  - - —  r  cos0  .  (A4) 

w  _  ,2  , 


37 


Substituting  expressions  (A2)-(A4)  into  Eq.  (Al)  yields  separate  monopole  and  dipole 
equations.  The  monopole  equation  for  A0  is  given  by 


1  9  3  „ 

-  r  T—  A  = 
r  or  or  o 


3A 

o  o 
c  3  r 


4  jt 

c  n 


where  jn  is  identified  as  the  net-current  density.  Solutions  to  this  equation  show  that  A0  relaxes 
on  the  L/R  time  discussed  in  Sec.  13-E. 


The  corresponding  dipole  equation  is 


3_  1  3_  _  4» 

3r  r  3r  1  c 


3  j  3A  ’ 

n  a  1 

*b  3r  c  3r  , 


where  A,  is  the  dipole  potential  given  by 


3A  21 

,  o  nr 

A  =  -  x  - - X - 

1  n  3r  n  a  c  a 

w  w 


2  .  ,  , 
x  —  i  (r  - 
n  rc  n 


Here  we  have  used  the  fact  that  the  monopole  magnetic  field  is  given  from  Eq.  (A5)  by 
B0  =  -bAJbr  =  2in(r)/rc,  where 


l  (r)  =  f  dr'  2»r'j 
n  Jn  ^ 


is  the  net  current  enclosed  within  radius  r  and  In  =  inf^).  Inserting  Eq.  (A7)  into  Eq.  (A6) 
produces,  after  some  manipulation. 


I  ,  ,21  -2oi  3x 

n  fr  )  _ n  _ n 

i  (r)  la  J  2  .  _  3t  ^  Xn 

n  w  ,  rc  8]  /3r. 

n 


where  we  have  taken  I„  to  be  constant  in  time. 

If  o  and  jn  have  the  same  Bennett  profile  of  radius  the  term 
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(A10) 
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rc^3j  /dr 
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do 


is  independent  of  r.  This  is  a  unique  property  of  Bennett  profiles,  and  it  justifies  treating  the 
dipole  field  as  a  simple  displacement  of  the  monopole  field,  plus  a  contribution  from  the  wall 
current. 

Equation  (A9)  shows  that  a  straw  reduces  the  relaxation  time  of  the  dipole  field  from  to 


r 


d 


r  i 

.  .21 

i  n  r 
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w  b 

do 


(All) 


for  Bennett  profiles.  The  dependence  on  r  is  weak  for  most  of  the  beam  if  a^,  £  3  a^  and  it  can 
be  "removed"  by  averaging  over  the  beam  to  define  an  effective  dipole  decay  time  given  for 
Bennett  profiles  by 


"V 
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3i, 

I  ,  21 
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1  -  n  r 

3r 

1  i  (r)  a 

0 

n  w 

'dol1  -  (Vaw'2  ln<1*aX>  ]  • 


(A12) 


For  flat-topped  profiles,  averaging  is  not  needed,  and  =  ( 1  -rb2/aw2),  where  rb  is  the  edge 
radius  given  by  Eq.  (30);  this  reduction  is  evident  in  the  analysis  of  Uhm  and  Lampe.39  Not 
surprisingly,  the  reduction  is  roughly  the  same  for  Bennett  profiles  when  expressed  in  terms  of 

rb- 

Let  us  now  add  orbit  phase  mixing.  Phase  mixing  occurs  when  the  local  pinch  force, 
proportional  to  i„(r)/r,  does  not  increase  linearly  with  r.  The  betatron  oscillation  frequency  of  the 
beam  electrons  then  varies  with  r,  attaining  its  maximum  value  k^0  on  the  beam  axis.  To  model 
this  variation  without  following  the  three-dimensional  trajectories  of  individual  electrons,  Lee35 


39 


represented  each  slice  of  the  beam  by  an  ensemble  of  components,  each  of  which  has  a  Bennett 
radial  profile  and  is  allowed  to  undergo  transverse  displacements  without  internal  distortion. 
However,  each  component  is  assigned  a  beam  electron  mass  which  is  weighted  by  a  variable 
V  *  1 .  so  that  a  component  oscillates  at  a  betatron  frequency  f»7  k^.  The  distribution  of  masses 
f(rj)  was  chosen  to  be  the  simplest  form  that  satisfies  certain  conservation  properties.  Overall, 
the  spread-mass  model  gives  good  results,  partly  because  the  hose  growth  rate  is  insensitive  to 
details  in  f(tj). 

Each  mass  component  separately  satisfies  Eq.  (31),  which  can  be  written  as 

a2  2 

~2  %  1  ”k0o  ( W  ■  lrb/aw 

02 

=  ”kfo  (W  -  (ab/aw 

where  rbks  =  abk^0  for  Bennett  profiles.  This  equation  applies  for  all  i?  between  0  and  1. 
Observe  that  we  did  not  apply  phase-mixing  to  the  wall  force  because  that  force,  unlike  the  self¬ 
pinch  force,  is  independent  of  position  and  is  therefore  the  same  for  all  particles  in  a  given  slice. 

To  derive  a  dispersion  relationship,  we  impose  solutions  of  the  form  x  -  e'kz',WT  on 
Eqs.  (A9)  and  (A  13).  That  reduces  Eq.  (A  13)  to 
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(A14) 


and  Eq.  (A9)  to 


(l-iur  )  x  , 
a  n 


(A15) 


where  rd  is  the  effective  dipole  decay  time  given  by  Eq.  (A  12).  We  then  integrate  Eq.  (A  14) 
over  f(r?)  and  combine  it  with  Eq.  (A15)  to  obtain 
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(A16) 
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where  x  =  k2/k^2  and  the  spread-mass  function  G  is  defined  by 


(A17) 


Equation  (A16)  is  the  spread-mass  dispersion  relationship  in  the  presence  of  wall  forces. 

The  functions  f(ij)  and  G(x)  depend  on  the  radial  structure  of  in(r).  A  straw  can  alter  in(r) 
by  restricting  the  maximum  orbit  radius  of  the  beam  electrons;  this  reduces  the  spread  in  orbit 
frequencies  and  forces  f(rj)  to  zero  at  small  17.  However,  the  net  effect  on  phase  mixing  and  G(x) 
is  not  severe,  provided  the  orbit  frequencies  vary  by  two  or  more  inside  the  straw.  This  is  easily 
met  for  Bennett  profiles  if  aw  ^  3  a^.  We  shall  therefore  impose  this  restriction  and  use  G(x)  as 
calculated  by  Lee  in  the  absence  of  wall  forces.  For  unbounded  Bennett  profiles,35 


G(x) 


6x  [j  -  X  +  X(l-X)  lnp^*)]  • 


(A18) 


If  growth  in  z  is  more  important  than  growth  in  r,  we  can  take  u  to  be  real.  The  real  part  of 
the  right-hand  side  of  Eq.  (A  16)  then  equals  zero,  yielding 


(a,  /a  )2 

D  W 


Re{G (x) } 
Re{G(X)/x) 


F(kr,k.)  . 


(A19) 


The  peak  growth  rate  k^*ax  occurs  where  3F/3kr  =  0.  This  allows  us  to  compute  Id-18*  as  a 
function  of  ab/aw.  In  Fig.  4,  we  plot  Kf^/k^  as  a  function  of  a^/a^,  using  the  Bennett  spread- 
mass  function  given  by  Eq.  (A  18).  Recall  that  this  function  improperly  includes  contributions 
from  outside  the  straw,  so  that  the  analysis  becomes  suspect  at  a^  <  3aj,.  For  comparison,  we 
also  show  the  flexible-rod  result  (dashed  curve).  The  good  agreement  between  the  two  curves 
indicates  that  phase-mixing  reduces  only  moderately. 

To  complete  the  analysis,  we  compute  the  peak  growth  rate  cJ[‘ax  in  r  for  real  k.  Recall  that 
cJ?8*  becomes  important  if  z/X^  »  Tp/rd.  For  real  k,  Eq.  (A16)  yields 

«±rd  =  [l  -  (ab/aw)2x"1]  Re{G(x>)  ,  <A20) 

where  x  =  (k/k^)2  is  real.  For  a  given  ab/aw,  the  peak  growth  is  easily  determined  and  is  plotted 
in  Fig.  8,  using  Eq.  (A12)  to  account  for  the  reduction  in  rd  due  to  the  straw.  The  growth  rate  is 
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everywhere  finite  because  phase-mixing  spreads  the  resonance  at  k  =  ks.  A  comparison  of 
Figs.  4  and  8  shows  that  a  straw  lowers  Id-**  more  than  cJJ13*,  for  aw  *  Sat,. 
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Rg.  1  -  Generic  welding  or  materials  processing  system  consisting  of:  a  pulsed  beam  generator,  a  focusing  and 
steering  section  (if  required);  an  intermediate  transport  section  (e.g.,  a  straw);  a  free  transport  zone;  various  vacuum/ 
gas  interfaces;  and  x-ray  shielding  surrounding  the  entire  unit. 


lb(kA) 


Fig.  Sa  -  Allowed  beam  current  and  energy  (shaded  region)  for  zt  =  10  cm, 
ao  =  0.1  cm,  aw  =  <»,  V2,  and  ghS  30 
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Fig.  5b  -  Allowed  parameters  for  z(  =  20  cm 
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Fig.  5c  -  Allowed  parameters  for  zt  =  30  cm 
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Fig.  6a  -  Allowed  parameters  at  zt=  30  cm  using  a  straw  of  radius  aw  =  3ao 
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Fig.  6b  -  Allowed  parameters  using  a  straw  filled  with  1  atm  of  He 
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-  Loss  of  beam  current  inside  a  straw 


54 


X,  Y  (cm) 


0  10  20  30 

z  (cm) 


Fig.  7d  -  Centroid  of  a  beam  injected  at  an  angle  of  dx/dz  =  10  mrad 
from  the  axis  of  the  straw 
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